Numerical Contractor Renormalization Method for Quantum Spin Models 
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We demonstrate the utility of the numerical Contractor Renormalization (CORE) method for 
quantum spin systems by studying one and two dimensional model cases. Our approach consists of 
two steps: (i) building an effective Hamiltonian with longer ranged interactions up to a certain cut-off 
using the CORE algorithm and (ii) solving this new model numerically on finite clusters by exact 
diagonalization and performing finite-size extrapolations to obtain results in the thermodynamic 
limit. This approach, giving complementary information to analytical treatments of the CORE 
Hamiltonian, can be used as a semi-quantitative numerical method. For ladder type geometries, 
we explicitely check the accuracy of the effective models by increasing the range of the effective 
interactions until reaching convergence. Our results in the perturbative regime and also away from 
it are in good agreement with previously established results. In two dimensions we consider the 
plaquette lattice and the kagome lattice as non-trivial test cases for the numerical CORE method. 
As it becomes more difficult to extend the range of the effective interactions in two dimensions, 
we propose diagnostic tools (such as the density matrix of the local building block) to ascertain 
the validity of the basis truncation. On the plaquette lattice we have an excellent description 
of the system in both the disordered and the ordered phases, thereby showing that the CORE 
method is able to resolve quantum phase transitions. On the kagome lattice we find that the 
previously proposed twofold degenerate S = 1/2 basis can account for a large number of phenomena 
of the spin 1/2 kagome system. For spin 3/2 however this basis does not seem to be sufficient 
anymore. In general we are able to simulate system sizes which correspond to an 8 x 8 lattice for 
the plaquette lattice or a 48-site kagome lattice, which are beyond the possibilities of a standard 
exact diagonalization approach. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.40.Cx 
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Low-dimensional quantum magnets are at the heart of 
current interest in strongly correlated electron systems. 
These systems are driven by strong correlations and large 
quantum fluctuations - especially when frustration comes 
into play - and can exhibit various unconventional phases 
and quantum phase transitions. 

One of the major difficulties in trying to understand 
these systems is that strong correlations often gener- 
ate highly non trivial low-energy physics. Not only the 
groundstate of such models is generally not known but 
also the low-energy degrees of freedom can not be iden- 
tified easily. Moreover, among the techniques available 
to investigate these systems, not many have the required 
level of generality to provide a systematic way to derive 
low-energy effective Hamiltonians. 

Recently the Contractor Renormalization (CORE) 
method has been introduced by Morningstar and 
Weinstem 1 . The key idea of the approach is to derive 
an effective Hamiltonian acting on a truncated local ba- 
sis set, so as to exactly reproduce the low energy spec- 
trum. In principle the method is exact in the low en- 
ergy subspace, but only at the expense of having a pri- 
ori long range interactions. The method becomes most 
useful when one can significantly truncate a local ba- 
sis set and still restrict oneself to short range effective 
interactions. This however depends on the system un- 
der consideration and has to be checked systematically. 
Since its inception the CORE method has been mostly 
used as an analytical method to study strongly correlated 



systemsSi^. Some first steps in using the CORE ap- 
proach and related ideas in a numerical framework have 
also been undertaken£*§i&. 

The purpose of the present paper is to explore the nu- 
merical CORE method as a complementary approach 
to more analytical CORE procedures, and to system- 
atically discuss its performance in a variety of low di- 
mensional quantum magnets, both frustrated and unfrus- 
trated. The approach consists basically of numerical ex- 
act diagonalizations of the effective Hamiltonians. In this 
way a large number of interesting quantities are accessi- 
ble, which otherwise would be hard to obtain. Further- 
more we discuss some criteria and tools useful to estimate 
the quality of the CORE approach. 

The outline of the paper is as follows: In the first sec- 
tion we will review the CORE algorithm in general and 
discuss some particularities in a numerical CORE ap- 
proach, both at the level of the calculation of the effective 
Hamiltonians and the subsequent simulations. 

In section [H] we move to the first applications on one- 
dimensional (ID) systems: the well known two-leg spin 
ladder and the 3-leg spin ladder with periodic bound- 
ary conditions in the transverse direction (3-leg torus). 
Both systems exhibit generically a finite spin gap and a 
finite magnetic correlation length. We will show that the 
numerical CORE method is able to get rather accurate 
estimates of the groundstate energy and the spin gap by 
successively increasing the range of the effective interac- 
tions. 
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In section lTTTl we discuss two-dimensional (2D) systems. 
As in 2D a long ranged cluster expansion of the interac- 
tions is difficult to achieve on small clusters, we will dis- 
cuss some techniques to analyze the quality of the basis 
truncation. We illustrate these issues on two model sys- 
tems, the plaquette lattice and the kagome lattice. The 
plaquette lattice is of particular interest as it exhibits 
a quantum phase transition from a disordered plaque- 
tte state to a long range ordered Neel antiferromagnet, 
which cannot be reached by a perturbative approach. We 
show that a range-two effective model captures many 
aspects of the physics over the whole range of param- 
eters. The kagome lattice on the other hand is a highly 
frustrated lattice built of corner-sharing triangles. For 
spin 1/2 it has been studied both numerically and ana- 
lytically and it is one best-known candidate systems for 
a spin liquid groundstate. A very peculiar property is 
the exponentially large number of low-energy singlets in 
the magnetic gap. We show that already a basic range 
two CORE approach is able to devise an effective model 
which exhibits the same exotic low-energy physics. For 
higher half-integer spin, i.e. S = 3/2, this simple effec- 
tive Hamiltonian breaks down; we analyze how to detect 
this, and discuss some ways to improve the results. 

In the last section we conclude and give some per- 
spectives. Finally three appendices are devoted to (i) 
the density matrix of local building block, (ii) the calcu- 
lation of observables by energy considerations and (iii) 
some general remarks on effective Hamiltonians coupling 
antifcrromagnetic half- integer spin triangles. 

I. CORE ALGORITHM 

The Contractor Renormalization (CORE) method has 
been proposed by Morningstar and Weinstein in the con- 
text of general Hamiltonian lattice models^. Later We- 
instein applied this method with success to various spin 
chain models^. For a review of the method we refer the 
reader to these original papersi^ and also to a peda- 
gogical article by Altman and Auerbach— which includes 
many details. Here, we summarize the basic steps before 
discussing some technical aspects which are relevant in 
our numerical approach. 

CORE Algorithm : 

• Choose a small cluster (e.g. rung, plaquette, tri- 
angle, etc) and diagonalize it. Keep M suitably 
chosen low-energy states. 

• Diagonalize the full Hamiltonian H on a connected 
graph consisting of N c clusters and obtain its low- 
energy states |n) with energies s n . 

• The eigenstates \n) are projected on the tensor 
product space of the states kept and Gram-Schmidt 
orthonormalized in order to get a basis \ip n ) °f di- 
mension M Na . As it may happen that some of the 
eigenstates have zero or very small projection, or 



vanish after the orthogonalization it might be nec- 
essary to obtain more than just M c exact eigen- 
states. 

• Next, the effective Hamiltonian for this graph is 
built as 

M Nc 

tlNc = ^2 £n\lpn)(lpn\- (1) 

71=1 

• The connected range- N c interactions h^ nn are de- 
termined by substracting the contributions of all 
connected subclusters. 

• Finally, the effective Hamiltonian is given by a clus- 
ter expansion as 

//'"•••' = ]T ^ + h ij + E h v<< + ■■■ ( 2 ) 

i (ij) (ijk) 

This effective Hamiltonian exactly reproduces the low- 
energy physics provided the expansion goes to infinity. 
However, if the interactions are short-range in the start- 
ing Hamiltonian, we can expect that these operators will 
become smaller and smaller, at least in certain situations. 
In the following, we will truncate at range r and verify 
the convergence in several cases. This convergence natu- 
rally depends on the number M of low-lying states that 
are kept on a basic block. In order to describe quanti- 
tatively how "good" these states are, we introduce the 
density matrix in section ITTT1 

When the number of blocks increases, a full diagonal- 
ization is not always easy and one is tempted to use 
a Lanczos algorithm in order to compute the low-lying 
eigenstates. In that case, one has to be very careful to 
resolve the correct degeneracies, which is known to be a 
difficult task in the Lanczos framework. In practice such 
degeneracies arise when the cluster to be diagonalized is 
highly symmetric. If the degeneracies are ignored, often 
a wrong effective Hamiltonian with broken SU(2) sym- 
metry is obtained. As a consequence we recommend to 
use specialized LAPACK routines whenever possible. 

In the present work we investigate mainly SU(2) in- 
variant Heisenberg models described by the usual Hamil- 
tonian 

H = JjjSj ■ Sj (3) 
(ij) 

where the exchange constants Jy will be limited to short- 
range distances in the following. As a consequence of the 
SU(2) symmetry, the total spin of all states is a good 
quantum number. This also has some effects when cal- 
culating the effective Hamiltonian. It is possible to have 
situations where a low energy state has a non-zero over- 
lap with the tensor product basis, but gets eliminated by 
the orthogonalization procedure because one has already 
exhausted all the states in one particular total spin sector 
by projecting states with lower energy. 
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(a) 




(b) 



FIG. 1: (a) 2-leg ladder. Basic block is a 2 x 2 plaquette. (b) 
3-leg torus with rung coupling J± and inter-rung coupling Jm . 



Once an effective Hamiltonian has been obtained, it is 
still a formidable task to determine its properties. Within 
the CORE method different routes have been taken in the 
past. In their pioneering papers Morningstar and Wein- 
stein have chosen to iteratively apply the CORE method 
on the preceding effective Hamiltonian in order to flow 
to a fixed point and then to analyze the fixed point. A 
different approach has been taken in Refs. |3j|4|: There 
the effective Hamiltonian after one or two iterations has 
been analyzed with mean-field like methods and interest- 
ing results have been obtained. Yet another approach - 
and the one we will pursue in this paper - consists of a 
single CORE step to obtain the effective Hamiltonian, 
followed by a numerical simulation thereof. This ap- 
proach has been explored in a few previous studie o 5 i 6 i 7 . 
The numerical technique we employ is the Exact Diago- 
nalization (ED) method based on the Lanczos algorithm. 
This technique has easily access to many observables and 
profits from the symmetries and conservation laws in the 
problem, i.e. total momentum and the total S z compo- 
nent. Using a parallelized program we can treat matrix 
problems of dimensions up to ~ 50 millions, however the 
matrices contain significantly more matrix elements than 
the ones of the microscopic Hamiltonian we start with. 



II. LADDER GEOMETRIES 

In this section, we describe results obtained on ladder 
systems with 2 and 3 legs respectively. 

We want to build an effective model that is valid from a 
perturbative regime to the isotropic case Jij = J = 1. We 
have chosen periodic boundary conditions (PBC) along 
the chains in order to improve the convergence to the 
thermodynamic limit. 



A. Two-leg Heisenberg ladder 




FIG. 2: Ground-state energy per site and spin gap of a 2 x L 
Heisenberg ladder using CORE method with various range r 
using PBC. For comparison, we plot the best known extrap- 
olations 10 with arrows. 



subspace is formed by the singlet ground-state (GS) and 
the lowest triplet state. 

Using the same CORE approach, Piekarewicz and 
Shepard have shown that quantitative results can be ob- 
tained within this restricted subspace^. Moreover, dy- 
namical quantities can also be computed in this frame- 
work*. 

Since we are dealing with a simple system, we can com- 
pute the effective models including rather long-range in- 
teractions (typically, to obtain range-4 interactions, we 
need to compute the low-lying states on a 2 x 8 lat- 
tice with Open Boundary Conditions which is feasible, 
although it requires a large numerical effort). It is de- 
sirable to compute long-range effective interactions since 
we wish to check how the truncation affect the physical 
results and how the convergence is reached. 

In a second step, for each of these effective models, we 
perform a standard Exact Diagonalization (ED) using 
the Lanczos algorithm on finite clusters up to N c = 12 
clusters (N = 48 sites for the original model). The GS 
energy and the spin gap are shown in Fig. [21 The use 
of PBC allows to reduce considerably finite-size effects 
since we have an exponential convergence as a function 
of inverse length. CORE results are in perfect agreement 
with known results and the successive approximations 
converge uniformly to the exact results. For instance, 
the relative errors of range-4 results are 10~ 4 for the GS 
energy and 10 -2 for the spin gap. This fast convergence 
is probably due to the rather short correlation length in 
an isotropic ladder (typically 3 to 4 lattice spacingsii). 



The 2-leg Heisenberg ladder has been intensively stud- 
ied and is known to exhibit a spin gap for all cou- 
plingslLiii. 

In order to apply our algorithm, we select a 2 x 2 pla- 
quette as the basic unit (see Fig. 0(a)). The truncated 



B. 3-leg Heisenberg torus 

As a second example of ladder geometry, we have stud- 
ied a 3-leg Heisenberg ladder with PBC along the rungs. 
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This property causes geometric frustration which leads to 
a finite spin-gap and finite dimcrization for all interchain 
coupling Jj , 12 ! 1 ? , contrary to the open boundary condi- 
tion case along the rungs, which is in the universality 
class of the Heisenberg chain. 

Perturbation theory : The simple perturbation the- 
ory is valid when the coupling along the rung (J±_) is 
much larger than between adjacent rungs (J\\)- In the 
following, we fix J±_ = 1 as the energy unit and denote 
a = J\\/J±. 

On a single rung, the low-energy states are the follow- 
ing degenerate states, defined as 
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where u> = exp(i27r/3). The indices L and R represent 
the momentum of the 3-site ring k y — 2n/3 and —2tt/3 
respectively. They define two chiral states which can be 
viewed as pseudo-spin states with operators f on each 
rung defined by 

r+\-R)=0 r+\-L) = \.R) 
t-\>R) = \>L) r-\-L) = 

t z \-R) = \\-R) t*\-L) = -\\-L) 

These states have in addition a physical spin 1/2 de- 
scribed by a. 

Applying the usual perturbation theory for the inter- 
rung coupling, one findsi2ii4 

Hpert = - T + fE^( 1 + ^ T i T i + 7 i~ 7 f)) ( 5 ) 
(ij) 

where N is the total number of sites. 

This effective Hamiltonian has been studied with 
DMRG and ED techniques and it exhibits a finite spin 
gap As = 0.28 J|| and a dimerization of the ground 
stately 

Here we want to use the CORE method to extend the 
perturbative Hamiltonian with an effective Hamiltonian 
in the same basis for any coupling. 

CORE approach : As a basic unit, we choose a single 3- 
site rung. The subspace consists of the same low-energy 
states as for the perturbative result (Eq. Q) which are 
4- fold degenerate (2 degenerate 5 = 1/2 states). We can 
apply our procedure to compute the effective interactions 
at various ranges, in order to be able to test the conver- 
gence of the method. 

First, we write down the range-2 contribution un- 
der the most general form which preserves both SU (2) 



(spin) symmetry and simultaneous translation or reflec- 
tion along all the rungs : 

H r=2 = Na + 5> Tf r/ + c {r+ T - + T r T +)) 

m 

+#i ■ aj( ai + &it?t? + c^t+tJ + rf t/)) (6) 

In the perturbative regime given in JSJ, the only non- 
vanishing coefficients are given by : — —1/4, a\ = a/3 
and Ci = 4a/3. 

The parameters of the effective Hamiltonian can be 
obtained and their dependence as a function of the inter- 
rung coupling a is shown in Fig. |3| We immediately 
see some deviations from the perturbative result since 
coefficients in panel (i) and (ii) are non-zero and become 
as important as the other terms in the isotropic limit. 
Surprisingly, we observe that ci follows its perturbative 
expression on the whole range of couplings whereas a\ 
deviates strongly as one goes to the isotropic case but 
does not change sign. 
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FIG. 3: CORE coefficients (see Eq.|SJ for two coupled trian- 
gles as a function of the inter-rung coupling a = J«/J±. The 
parameters were computed using range-2 CORE. The coeffi- 
cients in panef (iii) have been divided by their vaiues in the 
perturbative fimit. They therefore aif start at 1. 

In order to study how the physical properties evolve as 
a function of Jii/Jj_, we have computed the GS energy 
and the spin gap both for a small-coupling case and in the 
isotropic limit, up to range 5 in the effective interactions. 

Small interrung coupling : We have chosen J\\/J± = 
0.25 which corresponds to a case where perturbation the- 
ory should still apply. Using ED, we can solve the effec- 
tive models on finite lattices and on Fig. 0] we plot the 
scaling of the GS energy and of the spin gap as a function 
of the system length L. Even for this rather small value 
of J|i/J_l, our effective Hamiltonian can be considered as 
an improvement over the first order perturbation theory. 
Moreover, we observe a fast convergence with the range 
of interactions and already the range-3 approximation is 
almost indistinguishable from ED results. 
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FIG. 4: GS energy per site and spin gap for a 3 x L Heisenberg 
torus with J\\/J± = 0.25. Results are obtained using the 
CORE method at various range r. 



The estimated gap is 0.16J|| and correspond to a lower 
bound since ultimately the gap should converge exponen- 
tially to its thermodynamic value. Our value is consistent 
with the DMRG onei^ (~ 0.2Jii), and is already reduced 
compared to the strong coupling resuh— (Ag = 0.28Jy). 

Isotropic case : We apply the same procedure in the 
isotropic limit. As expected, the convergence with the 
range of interactions is much slower than in the pertur- 
bative regime. We show on Fig.0th.at indeed the ground 
state energy converges slowly and oscillates around the 
correct value. These oscillations come from the fact that, 
in order to compute range-r interactions, one has to study 
alternatively clusters with an even or odd number of 
sites. Since this system has a tendency to form dimers on 
nearest- neighbour bonds, it is better to compute clusters 
with an even number of sites. 

For the spin gap, we find accurate results even with 
limited range interactions. In particular, we find that 
frustration induces a finite spin gap ~ 0.11 Jy in that 
system. As in the previous case, this is a lower bound 
which is in perfect agreement with DMRG studj*i£. 

Moreover, we find that the singlet gap vanishes in the 
thermodynamic limit as 1/L 2 (data not shown). This sin- 
glet state at momentum n along the chains corresponds 
to the state built in the generalized Lieb-Schultz-Mattis 
argument. Here, the physical picture is a two-fold de- 
generate GS due to the appearance of spontaneous dimer- 
ization. 

Spinon dispersion relation for the spin tube : One of 
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FIG. 5: Same as Fig. |I]for the isotropic case Jii = J± = 1. 



the advantages of this method is to be able to get infor- 
mation on some quantum numbers (number of particles, 
magnetization, momentum. . . ) For example, the effec- 
tive Hamiltonian H e g still commutes with translations 
along the legs, with the total S 1 * * and t z so that we can 
work in a given momentum sector (k x ,k y ) with a fixed 
magnetization S l z ot . By computing the energy in each 
sector, we can compute the dispersion relation. 
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FIG. 6: Spinon dispersion relation (see text) as a function 
of longitudinal momentum (in units of 7r). We only plot the 
lowest branch corresponding to k y = ±27r/3. The odd lengths 
run from 5 to 13. The lines are guide to the eyes for an 
extrapolation on both sides of n/2. 

In order to try to identify if the fundamental excitation 
is a spinon, we compute the energy difference between 
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the lowest S = 1/2 state when the length is odd [L = 
2p + 1) and the extrapolated GS energy obtained from 
the data on systems with even length 2p and 2p + 2. The 
data are taken from CORE with range-4 approximation. 
On Fig. |HJ we plot this dispersion as a function of the 
longitudinal momentum, relative to the GS with L = 2p. 

We observe a dispersion compatible with a spinon- 
like dispersion, which is massive with a gap at tt/2 
~ 0.05 ~ As/2. This result is consistent with a pic- 
ture in which the triplet excitation As is made of two 
elementary spinons. With our precision, it seems that 
the spinons are not bound but we cannot exclude a small 
binding energy. 

We have on overall agreement with results obtained in 
the strong interchain coupling regime^. 

Therefore, with CORE method, we have both the ad- 
vantage of working in the reduced subspace and not being 
limited to the perturbative regime. Amazingly, we have 
observed that for a very small effort (solving a small clus- 
ter), the effective Hamiltonian gives much better results 
(often less than 1% on GS energies) than perturbation 
theory. It also gives an easier framework to systemati- 
cally improve the accuracy by including longer range in- 
teractions. 

For these models, the good convergence of CORE re- 
sults may be due to the fact that the GS in the isotropic 
limit is adiabatically connected to the perturbative one. 
In the following part we will therefore study 2D mod- 
els where a quantum phase transition occurs as one goes 
from the perturbative to the isotropic regime. 



III. TWO DIMENSIONAL SPIN MODELS 

In this section we would like to discuss the applica- 
tion of the numerical CORE method to two dimensional 
quantum spin systems. We will present spectra and ob- 
servables and also discuss a novel diagnostic tool - the 
density matrix of local objects - in order to justify the 
truncation of the local state set. 

One major problem in two dimension is the more elab- 
orate cluster expansion appearing in the CORE proce- 
dure. Especially our approach based on numerical diago- 
nalization of the resulting CORE Hamiltonian faces prob- 
lems once the CORE interaction clusters wrap around 
the boundary of the finite size clusters. We therefore 
try to keep the range of the interactions minimal, but 
we still demand a reasonable description of low energy 
properties of the system. We will therefore discuss some 
ways to detect under what circumstances the low-range 
approximations fail and why. 

As a first example we discuss the plaquette lattice 
[Fig. [7| (a)], which exhibits a quantum phase transi- 
tion from a gapped plaquette-singlet state with only 
short ranged order to a long range ordered antifcr- 
romagnetic state as a function of the interplaquette 
coupling 16 i 17 i 18 i 19 . We will show that the CORE method 
works particularly well for this model by presenting re- 
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FIG. 7: (a) The plaquette lattice. Full lines denote the pla- 
quette bonds J, dashed lines denote the inter-plaquette cou- 
pling J', (b) The trimerized kagome lattice. Full lines de- 
note the up-triangle J bonds, dashed lines denote the down- 
triangle coupling J'. The standard kagome lattice is recovered 
for J'/J = 1. 



suits for the excitation spectra and the order parame- 
ter. It is also a nice example of an application where 
the CORE method is able to correctly describe a quan- 
tum phase transition, thus going beyond an augmented 
perturbation scheme. 

The second test case is the highly frustrated kagome. 
lattice [Fig. 0(b)] with non-integer spin, which has been 
intensively studied for 5 = 1/2 during the last few 
j.(. s properties are still not entirely un- 
derstood, but some of the features are well accepted by 
now: There is no simple local order parameter detectable, 
neither spin order nor valence bond crystal order. There 
is probably a small spin gap present and most strikingly 
an exponentially growing number of low energy singlets 
emerges below the spin gap. We will discuss a convenient 
CORE basis truncation which has emerged from a per- 
turbative point of viewSiSiiiSi and consider an extension 
of this basis for higher non-integer spin. 



A. Plaquette lattice 

The CORE approach starts by choosing a suitable de- 
composition of the lattice and a subsequent local basis 
truncation. In the plaquette lattice the natural decom- 
position is directly given by the uncoupled plaquettes. 
Among the 16 states of an isolated plaquette we retain 
the lowest singlet [K = (0, 0)] and the lowest triplet 
[K = (TT,n)]. The standard argument for keeping these 
states relies on the fact that they are the lowest energy 
states in the spectrum of an isolated plaquette. 

As discussed in appendix [S] the density matrix of a 
plaquette in the fully interacting system gives clear in- 
dications whether the basis is suitably chosen. In Fig. [S] 
we show the evolution of the density matrix weights of 
the lowest singlet and triplet as a function of the inter- 
plaquette coupling. Even though the individual weights 
change significantly, the sum of both contributions re- 
mains above 90% for all J'/J < 1. We therefore consider 
this a suitable choice for a successful CORE application. 

A next control step consists in calculating the spec- 
trum of two coupled plaquettes, and one monitors which 




J'/J J'/J 



FIG. 8: Density matrix weights of the two most important 
states on a strong (J-bonds) plaquette as a function of J'/J. 
These results were obtained by ED with the original Hamil- 
tonian on a 4 x 4 cluster. 
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FIG. 9: Low energy spectrum of two coupled plaquettes. 
The states targeted by the CORE algorithm are indicated by 
arrows together with their SU (2) degeneracy. 



states are targeted by the CORE algorithm. We show 
this spectrum in Fig. ED along with the targeted states. 
We realize that the 16 states of our tensor product ba- 
sis cover almost all the low energy levels of the coupled 
system. There are only two triplets just below the 5 = 2 
multiplet which are missed. 

In a first application we calculate the spin gap for 
different system sizes and couplings J'/J. The results 
shown in Fig. 1101 indicate a reduction of the spin gap for 
increasing J'/J. We used a simple finite size extrapola- 
tion in l/N in order to assess the closing of the gap. The 
extrapolation levels off to a small value for J'/J > 0.6. 
The appearance of a small gap in this known gapless 
region is a feature already present in ED calculation of 
the original modeUS, and therefore not an artefact of our 



FIG. 10: Triplet Gap for effective system sizes between 20 
and 52 sites, as a function of the interplaquette coupling J'/J. 
For J'/J > 0.5 a simple extrapolation in l/N is also dis- 
played. These results compare very well with ED results on 
the original model fl9|l . 



method. It is rather obvious that the triplet gap is not 
a very accurate tool to detect the quantum phase tran- 
sition within our numerical approach. We will see later 
that order parameter susceptibilities are much more ac- 
curate. 

It is well known that the square lattice (J'/J = 1) is 
Neel ordered. One possibility to detect this order in ED 
is to calculate the so-called tower of excitation, i.e. the 
complete spectrum as a function of S(S + 1), S being 
the total spin of an energy level. In the case of standard 
collinear Neel order a prominent feature is an alignment 
of the lowest level for each S on a straight line, form- 
ing a so called "Quasi-Degenerate Joint States" (QDJS) 
ensemble^, which is clearly separated from the rest of 
the spectrum on a finite size sample. We have calculated 
the tower of states within the CORE approach (Fig. llljl . 
Due to the truncated Hilbert space we cannot expect to 
recover the entire spectrum. Surprisingly however the 
CORE tower of states successfully reproduces the gen- 
eral features observed in ED calculations of the same 
modeler (a) a set of QDJS with the correct degener- 
acy and quantum numbers (in the folded Brillouin zone); 
(b) a reduced number of magnon states at intermediate 
energies, both set of states rather well separated from the 
high energy part of the spectrum. While the QDJS seem 
not to be affected by the CORE decimation procedure, 
clearly some the magnon modes get eliminated by the 
basis truncation. 

In order to locate the quantum phase transition from 
the paramagnetic, gapped regime to the Neel ordered 
phase, a simple way to determine the onset of long range 
order is desireable. We chose to directly couple the order 
parameter to the Hamiltonian and to calculate general- 
ized susceptibilities by deriving the energy with respect 
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FIG. 11: Tower of states obtained with a range-2 CORE 
Hamiltonian on an effective N — 36 square lattice (9- 
site CORE cluster) in different reduced momentum sectors. 
The tower of states is clearly separated from the decimated 
magnons and the rest of the spectrum. 
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FIG. 12: Staggered moment per site as a function of the 
rescaled applied staggered field for the plaquette lattice and 
different values of J' /J. Circles denote the approximate cross- 
ing point of curves for different system sizes. We take the ex- 
istence of this crossing as a phenomenological indication for 
the presence of Neel LRO. In this way the phase transition is 
detected between 0.5 < J' c /J < 0.6, consistent with previous 
estimates. The arrows indicate curves for increasing system 
sizes. 



to the external coupling. This procedure is detailed in 
appendix [5] Its simplicity relies on the fact that only 
eigenvalue runs are necessary. Similar approaches have 
been used so far in ED and QMC calculations 29 30 . 

Our results in Fig. 1121 show the evolution of the stag- 
gered moment per site in a rescaled external staggered 
field for different inter-plaquette couplings J' and dif- 
ferent system sizes (up to 8 x 8 lattices). We note the 
appearance of an approximate crossing of the curves for 
different system sizes, once Neel LRO sets in. This ap- 
proximate crossing relies on the fact that the slope of 
rriL(hN) diverges with increasing N in the ordered phase 
in our case^ .. We then consider this crossing feature as 



an indication of the phase transition and obtain a value 
of the critical point J c /J = 0.55 ± 0.05. This estimate 
is in good agreement with previous studies using various 
methodsiiiSiiS. We have checked the present approach 
by performing the same steps on the two leg ladder dis- 
cussed in section 111 Al and there was no long range mag- 
netic order present, as expected. 



B. kagome systems with half-integer spins 



In the past 10 years many efforts have been devoted 
to understand the low energy physics of the kagome an- 
tiferromagnet (KAF) for spins l/2 20 i 21 i 22 i 23 i 24 . At the 
theoretical level, the main motivation comes from the 
fact that this model is the only known example of a two- 
dimensional Heisenberg spin liquid. Even though many 
questions remain open, some very exciting low-energy 
properties of this system have emerged. Let us sum- 
marize them briefly: (i) the GS is a singlet (S = 0) and 
has no magnetic order. Moreover no kind of more ex- 
otic ordering (dimer-dimer, chiral order, etc.) have been 
detected using unbiased methods; (ii) the first magnetic 
excitation is a triplet (S — 1) separated from the GS by 
a rather small gap of order J/20; (hi) more surprisingly 
the spectrum appears as a continuum of states in all spin 
sectors. In particular the spin gap is filled with an expo- 
nential number of singlet excitations: A/Qinglets ~ 1.15^; 
(iv) the singlet sector of the KAF can be very well repro- 
duced by a short-range resonating valence bond approach 
involving only nearest-neighbor dimers. 

From this point of view, the spin 1/2 KAF with its 
highly unconventional low-energy physics appears to be 
a very sharp test of the CORE method. The case of 
higher half- integer spins 5 = 3/2, 5/2... KAF is also of 
particular interest, since it is covered by approximative 
experimental realizations^!. Even if some properties of 
these experimental systems are reminiscent of the spin 
1/2 KAF theoretical support is still lacking for higher 
spins due to the increased complexity of these models. 

In this section we discuss in detail the range-two CORE 
Hamiltonians for spin 1/2 and 3/2 KAF considered as a 
set of elementary up-triangles with couplings J, coupled 
by down-triangles with couplings J' [see Fig. 0(b)]. The 
coupling ratio will be denoted by a = J' /J. Before go- 
ing any further into the derivation of the CORE effective 
Hamiltonian let us start with the conventional degenerate 
perturbation theory results. Note that in the perturba- 
tive regime these two approaches yield the same effective 
Hamiltonian. 

As described in Appendix [U] the most general two- 
triangle effective Hamiltonian involving only the two spin 
1/2 degrees of freedom on each triangle can be written 
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in the following form: 

H = Na (a) + ^ (6 (o;)T i .ey-f J -.ey (7) 

+ai(a)di.a 3 

+h (a)ai .<7j {n.eij ) (fj.eij ) 
+c 1 (a)a i .a j (Ti.e ij + fj.e i3 )). 

In the spirit of Mila's approach 2 ^ for spin 1/2 the first 
order perturbative Hamiltonian in a can easily be ex- 
tended to arbitrary half-integer spin 5: 



Kagome Lattice 



= -aa 3 (8) 

x (1 - 2(25 + \)fi.e a ) (1 - 2(25 + l)T 3 .e b ) 

and the coefficients of Q in the perturbative limit are 
given as a x (a) = f, h(a) = ^(25+ l) 2 , Cl (a) = 
'-(2S + 1) and a (a) = b (a) = 0. 



2a i 
9 



1. Choice of the CORE basis 

As discussed in the previous paragraph we keep the two 
degenerate 5 = 1/2 doublets on a triangle for the CORE 
basis. In analogy to the the plaquette lattice we calculate 
the density matrix of a single triangle embedded in a 12 
site kagome lattice for both spin 5 = 1/2 and 5 = 3/2, in 
order to get information on the quality of the truncated 
basis. The results displayed in Fig.ll3lshow two different 
behaviors: while the targeted states exhaust 95% for the 
5 = 1/2 case, they cover only w 55% in the 5 = 3/2 
case. This can be considered a first indication that the 
range-two approximation in this basis might break down 
for 5 > 1/2 half integer spin, while the approximation 
seems to work particularly well for 5 = 1/2, thereby 
providing independent support for the adequacy of the 
basis chosen in a related mean-field study 2 ^. 

We continue the analysis of the CORE basis by moni- 
toring the evolution of the spectra of two coupled trian- 
gles in the kagome geometry (c.f. Fig. I22|) as a function 
of the inter-triangle coupling J', as well as the states se- 
lected by the range-two CORE algorithm. The spectrum 
for the spin 5 = 1/2 case is shown in Fig. We note 
the presence of a clear gap between the 16 lowest states 
- correctly targeted by the CORE algorithm - and the 
higher lying bands. This can be considered an ideal case 
for the CORE method. Based on this and the results 
of the density matrix we expect the CORE range-two 
approximation to work quite well. 

We compare these encouraging results with the spec- 
trum for the spin 5 = 3/2 case displayed in Fig.^j] Here 
the situation is less convincing: very rapidly (J' / J > 
0.45) the low energy states mix with originally higher 
lying states and the CORE method continues to target 
two singlets which lie high up in energy when reaching 
J' / J — 1. We expect this to be a situation where the 
CORE method will probably not work correctly when 
restricted to range-two terms only. 
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FIG. 13: Density matrix weights of the different total spin 
states in a triangle of a 12 site kagome cluster with S = 1/2 
and 5 = 3/2 spins. These results are obtained for the homo- 
geneous case a = 1. 
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FIG. 14: Spectrum of two coupled triangles in the kagome 
geometry with S = 1/2 spins. The entire lowest band contain- 
ing 16 states is successfully targeted by the CORE algorithm. 



Based on the two-triangle spectra shown above we used 
the CORE algorithm to determine the coefficients of the 
general two-body Hamiltonian [Eqn. Q]. For an inde- 
pendent derivation, see Ref. [32]. The coefficients ob- 
tained this way are shown in Figs . 1 1 61 and 1 1 71 for 5=1/2 
and 5 = 3/2 respectively. In the limit a«l the coeffi- 
cients can be obtained from the perturbative Hamiltonian 
[Eqn. ijSJl]. There are two classes of coefficients in both 
cases: ao and bo are zero in the perturbative limit, i.e. 
they are at least second order in a. The second class of 
coefficients (ai, b%, c\) are linear in a. For improved visu- 
alisation we have divided all the coefficients in the second 
class by their perturbative values. In this way we observe 
in Fig. that coefficients b\ and c\ change barely with 
respect to their values in the perturbative limit. How- 
ever a\ has a significant subleading contribution, which 
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FIG. 15: Spectrum of two coupled triangles in the kagome 
geometry with 5 = 3/2 spins. The 16 states targeted by 
the CORE algorithm are indicated by the arrows and their 
degeneracies. 



FIG. 17: Coefficients of the CORE range-two Hamiltonian 
for two coupled 5 = 3/2 triangles. The coefficients in panel 
(iii) have been divided by their values in the perturbative 
limit. 
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FIG. 16: Coefficients of the CORE range-two Hamiltonian 
for two coupled S — 1/2 triangles. The coefficients in panel 
(iii) have been divided by their values in the perturbative 
limit. 



leads to a rather large reduction upon reaching the a = 1 
point. It does however not change sign. 

The situation for the 5 = 3/2 case in Fig. 1171 is differ- 
ent: while the coefficients b\ and c\ decrease somewhat, 
it is mainly a± which changes drastically as we increase a. 
Starting from 1 it rapidly goes through zero (a m 0.07) 
and levels off to roughly -7 times the value predicted by 
perturbation theory as one approaches a = 1. In this case 
it is rather obvious that this coefficient will dominate the 
effective Hamiltonian. We will discuss the implications 
of this behavior in the application to the 5 = 3/2 kagome 
magnet below. 

Let us note that the behavior of the ax coefficient is 
mainly due to a rather large second order correction in 
perturbation theory. Indeed we find good agreement 
with the values obtained in the perturbative approach 



of Ref. 



2. Simulations for S = 1/2 

After having studied the CORE basis and the effec- 
tive Hamiltonian at range two in some detail, we now 
proceed to the actual simulations of the resulting model. 
We perform the simulations for the standard kagome lat- 
tice, therefore a = 1. We will calculate several distinct 
physical properties, such as the tower of excitations, the 
evolution of the triplet gap as a function of system size 
and the scaling of the number of singlets in the gap. 
These quantities have been discussed in great detail in 
previous studies of the kagome S = 1/2 antiferromag- 

nr | .20.21. 22. 23.24 

First we calculate the tower of excitations for a kagome 
5=1/2 system on a 27 sites sample. The data is plot- 
ted in Fig. The structure of the spectrum follows 
the exact data of Ref. |2l| rather closely; i.e there is no 
QDJS ensemble visible, a large number of 5 = 1/2 states 
covering all momenta are found below the first 5 = 3/2 
excitations and the spectrum is roughly bounded from 
below by a straight line in 5(5 + 1). Note that the tower 
of states we obtain here is strikingly different from the 
one obtained in the Neel ordered square lattice case, see 
Fig. EH 

Next we calculate the spin gap using the range-two 
CORE Hamiltonian. Results for system sizes up to 48 
sites are shown in Fig. ^3 together with ED data where 
available. In comparison we note two observations: (a) 
the CORE range-two approximation seems to systemat- 
ically overestimate the gap, but captures correctly the 
sample to sample variations, (b) the gaps of the small- 
est samples (effective iV=12,15) deviate strongly from 
the exact data. We observed this to be a general fea- 
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FIG. 18: Tower of states obtained with a range-two CORE 
Hamiltonian on an effective N — 27 kagome lattice (9-site 
CORE cluster). There is a large number of low- lying states in 
each S sector. The symbols correspond to different momenta. 
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FIG. 20: Logarithm of the number of states within the mag- 
netic gap. Results obtained with the CORE range-two Hamil- 
tonian. For comparison exact data obtained in Refs. |2j 22] 
are shown. The dashed lines are linear fits to the exact diag- 
onalization data. 
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FIG. 19: Spin gap of the kagome S = 1/2 model on var- 
ious samples, obtained with the CORE method (range-two 
and three). Exact diagonalization result are also shown for 
comparison where available. 



ture of very small clusters in the CORE approach. In 
order to improve the agreement with the ED data we 
calculated the two CORE range-three terms containing 
a closed loop of triangles. The results obtained with this 
extended Hamiltonian are shown as well in Fig. 1191 These 
additional terms improve the gap data somewhat. We 
now find the CORE gaps to be mostly smaller than the 
exact ones. The precision of the CORE gap data is not 
accurate enough to make a reasonable prediction on the 
spin gap in the thermodynamic limit. However we think 
that the CORE data is compatible with a finite spin gap. 

Finally we determine the number of nonmagnetic ex- 
citations within the magnetic gap for a variety of system 
sizes up to 39 sites. Similar studies of this quantity in 
ED gave evidence for an exponentially increasing number 
of singlets in the gap2ii2£. We display our data in com- 



parison to the exact results in Fig. [201 While the precise 
numbers are not expected to be recovered, the general 
trend is well described with the CORE results. For both 
even and odd N samples we see an exponential increase 
of the number of these nonmagnetic states. In the case 
of N = 39 for example, we find 506 states below the first 
magnetic excitation. These results emphasize again the 
validity of the two doublet basis for the CORE approach 
on the kagome spin 1/2 system. 



3. Simulations for S — 3/2 

We have also simulated the CORE Hamiltonian ob- 
tained above for 5 = 3/2. While the energy per site 
is reproduced roughly, unfortunately the spectrum does 
not resemble an antiferromagnetic spin model, i.e. the 
groundstate is polarized in the spin variables. This fact 
is at odds with preliminary exact diagonalization data on 
the original 5 = 3/2 model^. We therefore did not pur- 
sue the CORE study with this choice of the basis states 
any further. Indeed, as suggested by the analysis of the 
density matrix and by the evolution of the spectrum of 
two coupled triangles, we consider this a breakdown ex- 
ample of a naive range-two CORE approximation. It is 
important to stress that the method indicates its failure 
in various quantities throughout the algorithm, therefore 
offering the possibility of detecting a possible breakdown. 

As a remedy in the present case we have extended the 
basis states to include all the 5 = 1/2 and 5 = 3/2 states 
on a triangle (i.e. keeping 20 out of 64 states). Compu- 
tations within this basis set are more demanding, but 
give a better agreement with the exact diagonalization 
results. At the present stage we cannot decide whether 
the breakdown of the 4 states CORE basis is related only 
the CORE method or whether it implies that the kagome 
5 = 1/2 and 5 = 3/2 systems do not belong to the same 
phase. 
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IV. CONCLUSIONS 

We have discussed extensively the use of a novel nu- 
merical technique - the so-called numerical Contractor 
Renormalization (CORE) method - in the context of low- 
dimensional quantum magnetism. This method consists 
of two steps: (i) building an effective Hamiltonian acting 
on the low-energy degrees of freedom of some elementary 
block; and (ii) studying this new model numerically on 
finite-size clusters, using a standard Exact Diagonaliza- 
tion or similar approach. 

Like in other real-space renormalization techniques the 
effective model usually contains longer range interac- 
tions. The numerical CORE procedure will be most effi- 
cient provided the effective interactions decay sufficiently 
fast. We discussed the validity of this assumption in sev- 
eral cases. 

For ladder type geometries, we explicitely checked the 
accuracy of the effective models by increasing the range 
of the effective interactions until reaching convergence. 
Both in the perturbative regime and in the isotropic case, 
our results on a 2-leg ladder and a 3-leg torus are in good 
agreement with previously established results. This rapid 
convergence might be due to the small correlation length 
that exists in these systems which both have a finite spin 
gap. 

In two dimensions, we have used the density matrix 
as a tool to check whether the restricted basis gives a 
good enough representation of the exact states. When 
this is the for the plaquette lattice or the S = 

1/2 kagome lattice, the lowest order range- two effective 
Hamiltonian gives semi-quantitative results, even away 
from any perturbative regime. For example we can suc- 
cessfully describe the plaquette lattice, starting from the 
decoupled plaquette limit through the quantum phase 
transition to the Neel ordered state at homogeneous cou- 
pling. Furthermore we can also reproduce many aspects 
of the exotic low-energy physics of the S = 1/2 kagome 
lattice. 

Therefore within the CORE method, we can have both 
the advantage of working in a strongly reduced subspace 
and not being limited to the perturbative regime in cer- 
tain cases. 

We thus believe that the numerical CORE method can 
be used systematically to explore possible ways of gener- 
ating low-energy effective Hamiltonians. An important 
field is for example the doped frustrated magnetic sys- 
tems, where it is not easy to decide which states are im- 
portant in a low-energy description, and therefore the 
density matrix might be a helpful tool. 

APPENDIX A: DENSITY MATRIX 

In this appendix we introduce the density matrix of a 
basic building block in a larger cluster of the fully inter- 
acting problem as a diagnostic tool to validate or inval- 
idate a particular choice of retained states on the basic 



building block in the CORE approach. 

In previous applications of the CORE method, the 
choice of the states kept relied mostly on the spectrum of 
an isolated building block. While this usually gives rea- 
sonable results it is not a clear a priori where to place 
the cut-off in the spectrum. 

The density matrix of a "system block" embedded in 
a larger "super block" forms a key concept in the Den- 
sity Matrix Renormalization Group (DMRG) algorithm 
invented by S.R. White in 1992 34 and is at the heart of 
its success. Based on this and related ideas^ we propose 
to monitor the density matrix of the basic building block 
embedded in a larger cluster and to retain these states 
exhausting a large fraction of the density matrix weight. 

Consider now a subsystem A embedded in a larger 
system B. Suppose that the overall system B is in state 
|\&) (e.g. the groundstate). We write the wavefunction 
as: 

|*> = X>*.&I°><8>I 6 >> (Al) 

a , 6 

where the sum index a runs over all states in A and index 
b over all states in B \ A. The density matrix p A of the 
subsystem A is then defined as 

= (A2) 

b 

The eigenvalues of p A denote the probability of finding 
a certain state a in A, given the overall system in state 
|*). 

Practically we calculate the groundstate of the fully 
interacting system on a medium size cluster by exact di- 
agonalization, and then obtain the density matrix of a 
basic building block, e.g. a four site plaquette. The den- 
sity matrix of a building block is a rather local object, so 
we expect that results on intermediate size clusters are 
already accurate on the percent level. The density ma- 
trix spectra shown in Figs . 151 and 1131 have been obtained 
in this way. In the models considered, a density matrix 
weight of the retained states of at least 90% yielded rea- 
sonable results within a range-two CORE approximation. 
It is possible to allow for a lower overall weight, at the ex- 
pense of increasing the range of the CORE interactions. 



APPENDIX B: OBSERVABLES IN THE 
NUMERICAL CORE METHOD 

The calculation of observables beyond simple en- 
ergy related quantities is not straightforward within the 
CORE method, as the observables need to be renormal- 
ized like the Hamiltonian in the first placed. 

A somewhat simpler approach for measurements of 
symmetry breaking order parameters consists in adding 
a small symmetry breaking field to the Hamiltonian (for 
a review see Ref. [30|). 



Let us denote O the extensive symmetry breaking op- 
erator, such that the order parameter is related to its GS 
average value m — 1/N(ipo\0\ipo). The occurence of a 
symmetry broken phase can be detected by adding this 
operator to the Hamiltonian : 

n{5) = h-so (bi) 

Since on a finite-size lattice the order parameter vanishes 
by symmetry for 6 — 0, the ground-state energy per site 
varies quadratically for small 8 : 

e(8) ~ e - ^Xo<5 2 , 

where xo is termed the corresponding generalized suscep- 
tibility. In that way the second derivative of the energy 
with respect to 8 at 8 = offers one possibility to detect 
a finite order parameter in the thermodynamic limi^i. 

We found that another possibility to conveniently track 
the presence of a finite order parameter is to measure 
directly m{8) in finite field 

m{8) = (*a|d|* s ) = de(<5)/d<5 

by the Hcllmann-Feynman theorem. When plotting m(8) 
as a function of the rescaled field N8 for various system 
sizes we observe an approximate crossing of the curves if 
there is a finite order parameter and no crossing in the 
absence of the order parameter. 

APPENDIX C: GAUGE INVARIANCE ON 
HALF-INTEGER SPINS KAGOME LIKE SYSTEMS 

In this appendix, we discuss half-integer spin Hamil- 
tonians with triangles as the unit cell. The ground state 
manifold of each unit cell is generated by the four degen- 
erate lowest states that can be built out of 3 half-integer 
S spins, namely the four Stot = 1/2 states. The idea of 
selecting these states as a starting point to describe the 
whole system low energy properties was originally intro- 
duced by Subrahmanyam for S = 1/2 25 on the kagome 
lattice and later used by Mila 23 . More recently it was 
reintroduced by Raghu et a$& for arbitrary half-integer 
S in the context of a chain of triangles. All these ap- 
proaches are pertubative and state that the triangle cou- 
plings J is much larger than the inter-triangle one J'. 

Here we would like to discuss some general properties 
of any effective Hamiltonian that can be derived either 
by perturbative methods or more sophisticated ones such 
as CORE. In particular, we would like to point out that 
a gauge invariance appears as a direct consequence of the 
state selection. 

To be more specific, let us label 1,2,3 the sites of the 
triangle (see Fig. l2~T)l . In order to build a total spin 1/2 
out of the three S, spins 2 and 3 couple into a 5 + £(1/2) 
with e = ±1. The coupling with the remaining site 1 
produces a spin 1/2 with chirality e = ±1. Note that 
this definition of chirality is equivalent to Eqs. for 
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FIG. 21: Definition of chirality e (see text for details). 




FIG. 22: The two-triangle problem, a is the coupling ratio 
J'/J. 

spin S = 1/2 up to a global unitary transform which is 
just a redefinition of the chirality quantification axis. 

In the following, the four selected spin-chirality states 
on a triangle i will be denoted as H^,^)). These 
states are the eigenstates of the z components of spin 
a and chirality f (both are spin 1/2 like operators) with 
T z \\Ei,Ui)) = (Ei/2)\\Ei,Ui)) and o- z \\e l ,v i )) = Vi\\si,Vi)). 

Let us now turn to the two-triangle problem. As it 
can be seen in Fig. 1221 the Hamiltonian is invariant un- 
der reflections with respect to the [xx') axis. Moreover, 
the reflection can be taken independently on each trian- 
gle. As a consequence, both chiralities (rf and r?) are 
conserved by the effective Hamiltonian and the r part is 
of the form 1 + a(rf + rj) + &T?r?. For any fixed value 
of (£j,£j), the total spin of the system is conserved and 
thus the spin part is SU(2) invariant. As a conclusion 
the most general two-triangle Hamiltonian allowed is of 
the form: 

[3t ■ 3j + c)(l + a(r? + r/) + 6r?r?). 

Gauge transformation: The form of the above Hamilto- 
nian is the consequence of the particular choice we made 
for labeling the sites of the triangle (see Fig. \'22li : site 1 
of triangle t\ couples to site 1 of triangle t%. Although 
this gauge was convenient for the calculation, in general 
this choice can not be made simultaneously on all cou- 
ples of triangles of the lattice. So, it is essential to derive 
the form of the Hamiltonian in a generic situation where 
site i — 1,2,3 of triangle t\ couples to site j = 1, 2, 3 of 
triangle t2- 

The unitary transformations involved in the redefini- 
tion of the coupling sequence (see Fig. l2*3*jl are covered by 
the 3j symbols of elementary quantum mechanics. The 
problem of 3 half-integer spins S coupled into a total spin 
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FIG. 23: Three ways of coupling the three spins S on a tri- 
angle into a total spin 1/2 state. Each construction is related 
to the two others by the 3j symbols (see text). 

1/2 occurs to be particularly simple and independent of 
S. The form of the general effective Hamiltonian then 
reads: 

= ' °3 + c ( a )) 

x [1 + a(a)(n ■ e a + fj ■ e b ) 
+b(a)(fi ■ e a ){fj ■ eft)], 

where e a , a — 1,2,3 are three coplanar normalized vec- 
tors in a 120° configuration (for example, e*i = (0,1), 
e 2 = (— v/3/2,-1/2) and e 3 = (y/3/2,-1/2) in the x - z 



plane) and a, b are the labels of the original spins cou- 
pling triangles i, and tj . 

The kagome lattice: In the particular geometry of 
the kagome lattice [see Fig. [7| (b)], each triangular unit 
cell is coupled to six other triangular cells, each cor- 
ner being coupled twice. As a consequence, for each 
cell the contribution involving only fi ■ e a factorizes into 
2fi ■ (e*i + e?2 + (3*3) = 0. The corresponding terms are 
then not relevant in the Hamiltonian and thus we de- 
note the most general two-triangle Hamiltonian for the 
kagome lattice as: 

H = Na {a) + ^ IM")^ ' ZijTj ' &ij 
(id) 

+ai(a)(Ti ■ ffj 

+bi{a)ai ■ aj(n ■ &ij)(Tj • e„) 
+ci(a)ai ■ c?j(fi ■ eij) + (fj ■ %)] 

which is the form used in the text. 
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